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I. INTRODUCTION 


A. OVERVIEW 


Ships can be initially designed for a variety of purposes or uses. 
But, as the times change, they may be no longer needed for this original 
use or may be considered for other types of operations which could 
require it to operate in a different environment than originally intended 
and designed for. The operating characteristics, in this new 


environment, may not be suitable for the intended use of the ship. 


The United States Air Force has been considering the use of T- 
AGOS class 3 and 4 ships as sea-based radar platforms. These ships 
were originally designed for surveillance and ASW operations. The 
operating characteristics for the original mission, and thus what the ship 
was designed for, are different to what would be required for this current 
consideration. To act as a sea-based radar platform, consideration 
would have to be given to improving the seakeeping qualities of the ship, 
in particular with regards to roll motion. Various methods exist for 


improving the roll of ships and are discussed in the next section. 


In this study, the effects of bilge keels on the ship were examined. 
In particular, various sizes of bilge keels for the full range of ship speed 


and sea states. A FORTRAN program, SHIPMO [Ref. 1], was utilized to 


perform the roll angle calculations for the varying conditions. The data 
obtained was utilized, in the form of various graphs and polar plots, to 
show the effects of bilge keels. Operability indices were found for each of 
the bilge keels under examination. If properly compared, these 
operability indices can provide a useful comparison of the effectiveness of 


bilge keels. 


B. STABILIZING METHODS 


Several methods exist for improving the seagoing stability of the 
ship and include, but not limited to, fin stabilizers, roll tanks, and bilge 
keels and may be used together or individually. The stability 
characteristics of the ship and the effects of employing these various 
methods can easily be seen by looking at the roll characteristics that are 
experienced by the ship. Various advantages and disadvantages exist for 
each of these methods and must be considered when determining a 


particular method to be used. [Ref. 2] 


The first type of stabilization technique that is available are active 
fin stabilizers. These produce a controlled roll moment where the phase 
and amplitude is such that it counteracts the external heel moment. The 
fins are most effective at higher speeds, generally greater than ten knots, 


since the force on the fin varies directly proportional to the speed of the 


ship. Below ten knots, the stabilizing moment available is not adequate 
to effectively counter the heel moment. In addition, the load on the fin or 
fin pivot may be such that the fin is unable to undergo full angular 
displacement at certain speeds. This could result in degraded 
performance at other additional speeds. Various parameters such as 
shapes, ship locations, and angles of attack can be varied to help obtain 
the maximum performance from them. The following are advantages and 


disadvantages of active fin stabilizers: 


Advantages: 
1. They offer the highest possible roll reduction with no reduction in 
static stability characteristics. They are the most effective of all single 


stabilizing devices. 


2. They are used in ships of different sizes. 


3. They inflict very small increases in ship resistance and have small 


auxiliary power requirements. 


Disadvantages: 


1. They are not effective at low speeds. 


2. They take up moderate machinery space, especially if they are 


retractable. This is desirable as they are less prone to damage. 


3. High initial cost due to the controlling equipment and machinery 


required. 


The next method of stabilization that exists are roll tanks. These 
tanks can be active or passive stabilizers which uses ship’s motion in 
such a way to cause water in the tanks to move in a direction as to 
oppose the ship’s motion. Two different tank configurations can be 
employed - free surface or U-tube. A major limitation to this method is 
that the tanks can only be tuned to one frequency. This is usually the 
natural frequency since this is most likely where the largest roll angles 
will occur. They are not as effective at other frequencies and can actually 
increase the roll angles attained at these other frequencies More 
sophisticated methods exist where the resistance in the duct can be 
adjusted for the frequency of the exciting waves which allows dampening 
to be achieved for all frequencies. The advantages and disadvantages for 


this method of roll stabilization are as follows: 


Advantages: 


1. High roll reduction rates of up to 70% are possible. 


2. They remain effective at low speeds. 


3. Vary little auxiliary power is required. 


4. Moderate initial cost and low maintenance required. 


5. Not vulnerable to damage since they are internal to the hull. 


Disadvantages: 


1. Moderate space is required in the hull. 


2. 1-4% reduction in deadweight capacity for constant displacement. 


3. Reduction in initial static stability due to free surface effects, 


especially for free surface tanks. 


The final method of stabilization that was considered was bilge 
keels. Bilge keels are the most widely used and simplest kind of roll 
stabilization in current use. They consist of a fin fixed to the hull at or 
near the bilge. The natural period of the roll of the ship is proportional to 
the radius of gyration of the ship. By attaching this fin, the radius of 
gyration of the ship is effectively increased. This results in an increased 


mass of water to roll with the ship and therefore an increase in the 


period of the roll. Under forced rolling conditions, such as in a seaway, 
the increased natural period that results from the bilge keel results in a 
roll amplitude reduction. The increased resistance to roll due to viscous- 
eddy effects supplied by the bilge keels plays an even larger role in roll 
reduction. Energy is dissipated by viscous flow from around the ship 
and this energy dissipation is increased substantially by bilge keel use. 


The advantages and disadvantages of this method are as follows: 


Advantages: 


1. Bilge keels are simple and easy to fit. 

2. They remain effective at relatively low speeds. 

3. Negligible reduction in ship’s deadweight capacity, no reduction in 
initial static stability, no auxiliary power requirements, and negligible | 
Space occupied by the hull. 

4, Low initial cost. 

Disadvantages: 


1. Since they are external to the hull, there is added resistance to ahead 


motion that must be overcome by the main engines. 


2. Comparatively to other methods, bilge keels offer smaller amounts of 


roll reduction. 


3. Bilge keels, carefully aligned to flow around the hull in calm waters to 
reduce forward motion resistance, can lead to added resistance during 


roll motion. 


4. Vulnerable to damage. 


All advantages and disadvantages of the various methods must be 
considered when determining which stabilization method is best suited 
for the application of interest. Some particular advantages or 
disadvantages may be of greater importance in the decision process than 
others. Active fin stabilizers, although the most effective, suffer from 
severe degradation at low speeds. Passive roll tanks have comparatively 
large response times and quick changes in roll motion can disrupt the 
operation of the roll tanks. Active roll tanks require large initial costs 
and auxiliary power with minimal added advantage. Due to the 
simplicity and relative effectiveness at low speeds, bilge keels were the 


method of choice. 





Il. EVALUATION OF RESPONSE 


A. REGULAR WAVE MOTION 


An important consideration for any floating body is the effects of 
waves on the free surface, especially surface waves with a period of a few 
seconds. Other waves, such as subsurface waves, usually have much 
longer periods and have negligible effect compared to surface waves. 
When analyzing these types of waves, some assumptions must be made. 
These include that the fluid is ideal and body motions are sufficiently 
small to linearize. From this, the appropriate fluid mechanics tools can 
be utilized to describe sea waves and ship motions based on potential 
flow theory. [Ref. 3] 

Potential, or ideal, flow theory is based on a couple of assumptions. 
The first assumption is that mass is preserved. If a control volume is 
placed around the object of interest and the fluid around it, then the 
mass which enters the volume must either accumulate in the control 
volume or leave. Utilizing the divergence theorem, conservation of mass, 
and assuming that density remains constant, the final form of the 


continuity equation inside the control volume becomes: 


—+—+—=0 (1) 


The next assumption is that flow is irrotational. This results in the 


property that circulation around any closed curve is Zero or: 


[U-aF=0 (2) 


From equation (2) and the definition of a velocity potential, the continuity 


equation can be written as: 


2 a cate, (3) 


aor’ oy? Oe? ve 
which is better known as Laplace’s Equation. 
Once the velocities and ¢have been solved for and assuming 
unsteady, irrotational flows, pressure can be found with the help of 


Bernoulli's equation and can be computed from: 


ee ee 
p=—p--—{ PV: Vo— paz (4) 


The plane progressive wave system is the simplest free surface 
wave formation scheme. It is two dimensional, sinusoidal in time, and 
propagates with a phase velocity such that an observer moving with the 
same velocity will make the wave appear stationary. Adopting a 
Cartesian coordinate system, the free surface elevation can be expressed 
in the following form: 

n(x,t) = Acos(kx — at) (5) 
where: A-wave amplitude 


k-wave number 
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Ship motions induced by regular wave motion are of most concern. 
In the simplest case, the waves incident upon the body may be assumed 
as plane progressive waves of small amplitude and sinusoidal time 
dependence. Ship motion is also assumed sufficiently small for linear 
theory to hold. Waves incident to a stationary or moving body will cause 
the body to move with six degrees of freedom or to heave, pitch, roll, yaw, 
surge, and sway. In general, the body motion in any of the six degrees of 


freedom can be expressed as: 


S -] 
nj =A>|D,| F (6) 
j=l 
where: A-wave amplitude 


F-exciting force per unit wave amplitude 
D-matrix of equations of motion 
A new fundamental quantity can also be introduced, known as the 


Response Amplitude Operator (RAO), as follows: 
DOW Os” (7) 
J 9 ? A 


This corresponds to the complex amplitude of body motion in the jth 
mode in response to an incident wave of unit amplitude, frequency, and 
direction, and is generally known as the transfer function, which can be 
calculated once the added mass, damping, exciting, and hydrostatic 


forces are known. 


bl 


The responses derived from the above motions are important for 
evaluating the seakeeping performance of the ship. Computation of the 
RAO’s is simple once the added mass and damping coefficients are 
known. But, this is the most difficult part. A technique for determining 
these values utilizes the Strip theory. In this method, the ship is 
assumed as a slender body. In other words, the beam and draft are 
much smaller than the length of the ship. The cross-section of the ship 
are divided into “strips”, or as for Frank’s mete [Ref 4]- a series of 
straight line segments integrated along the entire length of the ship so as 
to determine the overall effect. By utilizing the strip theory method, the 
ship’s motions in regular waves can be determined. Through a 
combination of the above equations, the resulting linearized 
simultaneous equations that must be solved for to determine these 


motions are: 


Yi 4-2 (Mp + Ay )+io,B, + Cy ra =F'+F?  j=1..6 (8) 
where: M - mass matrix 

A-added mass matrix 

B-dampening matrix 

C-hydrostatic restoring force matrix 


F' -Froude Krylov exciting force in the jth mode of motion 


F° -Diffraction exciting force in the jth mode of motion 
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¢, complex amplitude of motion in the kth direction 

B. ROLL DAMPING 

Roll damping can be accomplished by a variety of methods as 
discussed above. The roll erties hydrodynamic moment used in 
equations of motion tends to be nonlinear in nature with respect to the 
roll angle. To allow roll damping to be calculated, as in SHIPMO [Ref 5], 
it must first be linearized. Utilizing the method described by Himeno 
(1981) [Ref 6], the roll damping coefficient can be assumed in regular 


Waves aS: 
8 
By, (a =—b, + ee (9) 


where: B - damping matrix/coefficients 
¢ - Complex amplitude of motion 


@ - wave frequency 


The damping coefficients are a function of frequency, the point 
about which roll occurs, and the forward speed. In regular wave motion 
with constant forward speed, the damping coefficient values are constant 
about a fixed point. When irregular or random waves are encountered, 


the following equation by Himeno applies: 
8 
B,,(C,) = B,+ — OB: (10) 


where: oO - variance in angular roll velocity 
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The damping coefficients are composed of the following 


components: 
B, = B,+B,+B, 11) 
5, = Bo Ba: 

where: B, - frictional damping component 


By - wave generation damping component 

B, - hull lift damping component 

B_ - hull eddy damping component 

Bz, ~ bilge keel damping component 

The above components are based mainly on experimental analysis 

which are measured about a fixed point (roll center) and are usually 
limited in their range of applicability. If the beam to draft ratios are 
much more than about 2.5, the eddy dampening by Himeno tends to 
overpredict the eddy component. In addition, determining the point 
about which the ship is rolling, or roll eree can ral to be difficult. In 


particular, a ship in a seaway has only an instantaneous roll center and 


can be calculated as follows: 


= ig —a,)—tan(@,t+a,)sin(a, - a,)} (12) 
4 
where: a,- sway phase angle 


a, - roll phase angle 


14 


C. MOTIONS IN A SEAWAY 


Wave patterns in an open sea are ever changing with respect to 
time and space, in a manner that appears to defy analysis whether it is 
linear or second order Stokes. [Ref 7 ] Ambient waves on the surface of 
the sea are dispersive as well as random in nature where random refers 
to the character of the wave height distribution. The distribution of 
sinusoidal waves are such that they have continuously distributed 
amplitudes and phases such that the summation of the variation of wave 
height with time is not systematic in any respect, but random. The 
generating mechanism is, predominantly, the effect of wind in the 
atmosphere upon the water surface. Spectral Density, S(@), provides 
useful data that has been removed from a random wave record, h(t). The 
random wave record can be processed such that a S({@) vs. wave 
frequency, @, curve can be generated. The spectral density is obtained 
from a wave height record taken during a period of time when sea 
conditions are considered unchanging (stationary) for a certain sea state. 
The spectral energy density, S(@,@), (or directional energy spectrum) can 
be integrated over all wave directions, 0, to give the frequency spectrum 


as given below: 


S(@) = | S(@,a)de (13) 
0 
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For open seas, it is appropriate to assume that waves are unidirectional which 
allows the energy spectrum to be proportional to a delta function in 8 where the wave 
crests can be considered parallel and the fluid motion to be two dimensional. Wave 
spectra of this form is considered long crested and is sufficient to describe the wave 


environment. 


In naval architecture and ocean engineering, it is usually appropriate to assume 
waves to be long crested. With this simplification, existing information for the energy 
spectrum, based on theory and full scale observations, can be used for the energy 
spectrum. Usually, we are more concerned with larger waves. The most common 
parameter that takes this into account is the significant wave height, H,,,, which is 
defined as the average of the highest one third of all waves and can be expressed as 
follows: 

Fi a0 (14) 
Where m, is the integral of S(@ ) over all frequencies. More generally, the moments, 


m,, of the spectrum can defined by: 


m, = |o'S(@)do,i = 0,1,2,... (15) 
0 


Also, the average frequency of the spectrum, defined as the expected number of zero 


upcrossings of the body that is considered per unit time (i.e. the number of times the 
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wave amplitude passes through zero with a positive slope), can be 


determined as follows: 





o,=(—)" (16) 
mM, 
The average period between zero upcrossings is: 
le = 2n = Qn mM, (i 7) 
QO. mM, 


H1/3=1to8m 


S(w) (m2sec) 





w (1/sec) 
Figure 1. Typical Pierson-Moskowitz wave spectrum 
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A good model for fully developed seas is the Pierson-Moskowitz 
spectrum. Figure ]1 shows a typical Pierson-Moskowitz spectrum. This 
spectrum is based on significant wave height where it is assumed that 
the wave spectra has reached a steady-state equilibrium and is 
independent of duration and fetch. In other words, the waves are 
considered fully developed. It is recognized as an asymptotic form 
reached only after an extended period of wind with no contamination 


from an underlying swell and can be given as follows: 
S(w) = 0.00812? /@° exp(-0.032(g / H,o*) | (18) 


where: g - acceleration of gravity 
H, - significant wave height 


@ - wave frequency 


Once the frequency spectrum is found, another een can be 
determined which will allow the spectral description of waves to be 
generalized into regular harmonic waves. As previously discussed, the 
Response Amplitude Operator, RAO, is a nondimensional quantity that 
can be used to relate the amount of heave amplitude per unit wave 
length. If the sea waves are described by a random distribution and if 
body response to each wave is defined by a RAO, the body response will 


be: 
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n,(t) = 8 |[Z,@,0d4@0) (19) 
where: KR - any body response 

Z - response amplitude operator 

@ - wave frequency 


8 - wave direction 


The body response is also a random variable. The seaway spectrum may be 


related to any body response by the following equation: 
Sp) =|Z,@)2S@) (20) 


where: S(@ ) - spectrum of the seaway 


The above RAO is not only valid for regular wave motions but also in a 
spectrum of random waves. This allows for the assumption that a vessel with favorable 
characteristics in regular waves will also have the same favorable characteristics in 
irregular waves. From equation 20, various ship’s motions can be calculated from a 


given sea spectrum, such as the Pierson-Moskowitz discussed previously. 


i 


D. CALCULATION PROGRAM 


The computer program SHIPMO [Ref 1] was utilized to make the 
various calculations, more specifically roll calculations, needed for this 
study. The computer program calculates a variety of motion responses 
and is able to predict ship motions in six degrees of freedom. Regular 
wave motions are calculated utilizing a modified strip theory of Salvesen, 
Tuck, and Faltinsen (1970) [Ref 8] which is based on a slender body in 
an ideal fluid with corrections made for eae and other effects . Two 
dimensional properties are calculated using Frank’s method (1967) [Ref 
9] with a lid applied to remove any irregular frequencies which can cause 
calculation difficulties if not applied. Frank’s method is computationally 
fast and any ship’s cross-section can be approximated with as much 
accuracy as desired. The wave spectrum can either be user inputted or 
chosen from one of six methods available in the program. All 
calculations, for this study, were based on the classical Pierson- 
Moskowitz spectrum as described in the previous section. 

The program reads in the initial ship data and calculates the 
potential for the four modes of motion - pitch, roll, yaw, and sway, over 
the specified frequency range at each designated cross section utilizing 
Frank’s method. The calculations are based at the ships’ amidships 
centerline position unless otherwise specified. The program allows for a 


maximum of 21 stations to be selected along the length of the ship, with 
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a maximum of 15 input points for each station. These are required to provide the hull 
shape data necessary for the various calculations. The wave spectral ordinates are 
calculated by the Pierson-Moskowitz method as discussed above. The spectral 
ordinates are also used to determine the response amplitude operator as was also 
previously discussed. The motion amplitude is calculated in all six degrees of freedom 
with the horizontal and vertical plane motions solved separately. Roll dampening is 
calculated by the Himeno method with eddy dampening for box barges with sharp 
corners being calculated by either Ikeda’s (1977) [Ref 10] or Yamashita and Katagiri’s 
method (1980). [Ref 11] For bilge keel damping, the Himeno method is also utilized. 
From this, the desired ships motions can be calculated over ship’s speed, heading 


angle, and wave frequency variations. 


The program required a data initialization file, SHIPMO.IN, and is included in 
Appendix A. The file inputs a variety of constants required for the proper 
computations to be made. Bilge keel width, depth, and location are inputted from this 
file. Station positions make up the majority of the file which include the station 
position relative to the center of the ship. For each station, hull depth and width 
positions relative to the center and water lines, are provided with the first input point 


located at the keel of the ship. 
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E. SHIP AND ENVIRONMENTAL CONSIDERATIONS 


The analysis performed in this study were based on the Class 4 T-AGOS ship. 
Table 1 provides some of the general characteristics of the ship that were necessary for 


the calculations. 


Length 224 feet 


Location of Centerline 97 feet aft of bow 


Design Draft ISP Teet 





Maximum Ship Speed 11.3 knots 


Table 1. General Characteristics T-AGOS Class 4 Ships 


SHIPMO, as described in the previous section, was utilized to determine the 
roll characteristics of the ship based on different bilge keel characteristics. Hull width 
and depths for various locations along the length of the ship were required for the 
program’s calculations. These were obtained from the body plans for the T-AGOS 4 
class ship. [Ref 12] A scaled version of the one actually used may be viewed in 
Appendix B. The actual data inputted for each of the reference points may be found in 


the data input file SHIPMO.IN in Appendix A. 


Roll calculations were made through the entire speed and sea state spectrum. 


Speeds of the ship were incremented from zero to 12 knots in two knot speed 
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increments. Sea states are based on and provide an indicator of the 
heights of the waves, usually the significant wave height, being 
encountered by a vessel. They range from 1, relatively calm seas, to 8 
which are usually only encountered during strong storms such as 
hurricanes. Table 2 provides the significant wave height for each of the 


sea states considered. 


Sea State Significant Wave Height (feet) 






— 
61.6 


Table 2. Sea State vs. Significant Wave Height 
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lil. ROLL RESPONSE RESULTS 


A. ROLL CALCULATIONS 


As discussed, SHIPMO was utilized to make the roll calculations 
necessary for this study and a sample input file is included in Appendix 
A. The ship’s roll characteristics were analyzed based on several bilge 
keel sizes including no bilge keel. Appendix C provides a sample of the 
program that was utilized to analyze the data obtained and print out the 
results. Table 3 lists the different bilge keel sizes that were analyzed. 
The maximum bilge keel width was based upon two-thirds the length of 
the ship and subsequently divided in half for each smaller size. A 
maximum width of three feet was set and reduced in size in one foot 
increments. Each bilge keel was analyzed such that it was assumed 


centered along the centerline of the ship. 
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Length x Width 


(ft x ft) 





Table 3. Summary of bilge keel sizes analyzed 


To ensure that a good representative data set was obtained for each bilge keel 
configuration, enough data points were required. For each ship’s speed analyzed, the 
roll angle was found for each sea state at 22.5 degree increments of wave direction 


around a 360 degree circle centered on the ship. 


Once this data was obtained, a curve of significant roll angles placed on the ship 


for each sea state and a given speed and bilge keel configuration, could be found 


utilizing the spline curve fitting Ofunction of MATLAB [Ref 13] for each degree of 
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wave motion from O to 360 degrees. Enough data was required to 
ensure a smooth fit of the curves. To illustrate this point, a graph of the 
roll angles for no bilge keel and a ship’s speed of 2 knots may be seen in 
Figure 2. Similar graphs for all the other ship’s speeds and bilge keel 


configurations could also be produced if needed. 


Roll angle vs. wave angle (2 kts., sea states (ss) 1 to 8) 
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Figure 2. Roll Angle vs. wave angle (no bilge keel and 2 knots ship speed) 
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From the data obtained as just discussed, a splined curve of roll 
angles for all sea states for each degree of wave motion was found. From 
this, the maximum allowed sea states to produce a roll angle of 2, 4, or 6 
degrees were found. A polar plot of this data for each bilge keel 
configuration and respective speeds are shown in figures 3 to 72. In 
addition, polar plots combining all ship speeds for a certain bilge keel 
configuration were generated and can be observed in figures 73 to 102. 
The next section discusses the effectiveness of bilge keels. Operability 


indices were also determined and plotted to help show this. 


By examining the figures, several observations become evident. As 
was expected, the maximum allowed sea state increased as the 
maximum allowed roll angle was increased and as ship’s speed was 
increased, the maximum allowed sea state correspondingly decreased. 
Beam seas seemed to have the least effect on the roll characteristics of 
the ship and roll was assumed zero for sea directions relating to the bow 
or stern of the ship. The largest effect on roll angle came from quartering 
seas with quartering seas aft the beam having the largest effect. Thus, 
aft quartering seas provided the limiting position for roll angle limits. As 
can be seen, significant roll reduction can be expected by the use of bilge 


keels. 
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Figure 3. Sea state vs. wave angle for no bilge keel - 0 kts. 
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Figure 4. Sea state vs. wave angle for no bilge keel - 2 kts. 
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Figure 5. Sea state vs. wave angle for no bilge keel - 4 kts. 
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Figure 6. Sea state vs. wave angle for no bilge keel - 6 Kts. 
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Figure 7. Sea state vs. wave angle for no bilge Keel - 8 kts. 
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Figure 8. Sea state vs. wave angle for no bilge Keel - 10 kts. 
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Figure 9. Sea state vs. wave angle for no bilge keel - 12 kts. 
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Figure 10. Sea state vs. wave angle for 150x3 bilge keel - O kts. 
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Figure 11. Sea state vs. wave angle for 150x3 bilge keel - 2 kts. 
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Figure 12. Sea state vs. wave angle for 150x3 bilge keel - 4 kts. 
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Figure 13. Sea state vs. wave angle for 150x3 bilge keel - 6 kts. 
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Figure 14. Sea state vs. wave angle for 150x3 bilge keel - 8 kts. 
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Maximum sea state vs. wave angle (for 10 kt. and 2,4,6 degree roll angle) 
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Figure 15. Sea state vs. wave angle for 150x3 bilge keel - 10 kts. 
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Figure 16. Sea state vs. wave angle for 150x3 bilge keel - 12 kts. 
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Figure 17. Sea state vs. wave angle for 150x2 bilge keel - O kts. 
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Figure 18. Sea state vs. wave angle for 150x2 bilge keel - 2 kts. 
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Figure 19. Sea state vs. wave angle for 150x2 bilge keel - 4 kts. 
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Figure 20. Sea state vs. wave angle for 150x2 bilge keel - 6 kts. 
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Maximum sea state vs. wave angle (for 8 kt. and 2,4,6 degree roll angle) 
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Figure 21. Sea state vs. wave angle for 150x2 bilge keel - 8 kts. 
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Figure 22. Sea state vs. wave angle for 150x2 bilge keel - 10 kts. 
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wave angle for 150x2 bilge keel - 12 kts. 
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Figure 24. Sea state vs. wave 
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wave angle for 150x1 bilge keel - 2 kts. 


Figure 25. Sea state vs. 
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Figure 26. Sea state vs. 
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wave angle for 150x1 bilge keel - 6 kts. 
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Figure 27. 
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vs. wave angle for 150x1 bilge keel - 8 kts. 


Figure 28. Sea state 
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vs. wave angle for 150x1 bilge keel - 10 kts. 
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Figure 31. Sea state vs. wave angle for 75x3 bilge keel - O kts. 


Maximum sea state vs. wave angle (for 2 kt. and 2,4,6 degree roll angle) 
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Figure 32. Sea state vs. wave angle for 75x3 bilge keel - 2 kts. 
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Figure 33. Sea state vs. wave angle for 75x3 bilge keel - 4 kts. 
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Figure 34. Sea state vs. wave angle for 75x3 bilge keel - 6 kts. 
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Figure 35. Sea state vs. wave angle for 75x3 bilge keel - 8 kts. 
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Figure 36. Sea state vs. wave angle for 75x3 bilge keel - 10 kts. 
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Figure 37. Sea state vs. wave angle for 75x3 bilge keel - 12 kts. 
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Figure 38. Sea state vs. wave angle for 75x2 bilge keel - O kts. 
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Figure 39. Sea state vs. wave angle for 75x2 bilge keel - 2 kts. 
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Figure 40. Sea state vs. wave angle for 75x2 bilge keel - 4 kts. 
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Figure 41. Sea state vs. wave angle for 75x2 bilge keel - 6 kts. 
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Figure 42. Sea state vs. wave angle for 75x2 bilge keel - 8 kts. 
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Figure 43. Sea state vs. wave angle for 75x2 bilge keel - 10 kts. 
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Figure 44. Sea state vs. wave angle for 75x2 bilge keel - 12 kts. 
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Figure 45. Sea state vs. wave angle for 75x1 bilge keel - O kts. 
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Figure 46. Sea state vs. wave angle for 75x1 bilge keel - 2 kts. 
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Figure 47. Sea state vs. wave angle for 75x1 bilge keel - 4 kts. 
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Figure 48. Sea state vs. wave angle for 75x1 bilge keel - 6 kts. 
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Figure 49. Sea state vs. wave angle for 75x1 bilge keel - 8 kts. 
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Figure 50. Sea state vs. wave angle for 75x1 bilge keel - 10 kts. 
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Figure 51. Sea state vs. wave angle for 75x1 bilge keel - 12 kts. 
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Figure 52. Sea state vs. wave angle for 37.5x3 bilge keel - O kts. 
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Figure 53. Sea state vs. wave angle for 37.5x3 bilge keel - 2 kts. 
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Figure 54. Sea state vs. wave angle for 37.5x3 bilge keel - 4 kts. 
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Figure 55. Sea state vs. wave angle for 37.5x3 bilge keel - 6 kts. 


Il angle) 


ngle (for 8 kt. and 2,4,6 degree ro 





Figure 56. Sea state vs. wave angle for 37.5x3 bilge keel - 8 kts. 
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Figure 57. Sea state vs. wave angle for 37.5x3 bilge keel - 10 kts. 
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Figure 58. Sea state vs. wave angle for 37.5x3 bilge keel - 12 kts. 
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Figure 59. Sea state vs. 


wave angle for 37.5x2 bilge keel - 2 kts. 
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Figure 60. Sea state vs. 
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vs. wave angle for 37.5x2 bilge keel - 4 kts. 
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Figure 62. Sea state vs. wave angle for 37.5x2 bilge keel - 6 kts. 
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Figure 63. Sea state vs. wave angle for 37.5x2 bilge keel - 8 kts. 
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Figure 64. Sea state vs. wave angle for 37.5x2 bilge keel - 10 kts. 
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Figure 65. Sea state vs. wave angle for 37.5x2 bilge keel - 12 kts. 
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Figure 66. Sea state vs. wave angle for 37.5x1 bilge keel - 0 kts. 
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Figure 67. Sea state vs. wave angle for 37.5x1 bilge keel - 2 kts. 


Figure 68. Sea state vs. wave angle for 37.5x1 bilge keel - 4 kts. 





ngle (for 6 kt. and 2,4,6 degree roll angle) 


Maximum sea state vs. wave a 


90, 





Figure 69. Sea state vs. wave angle for 37.5x1 bilge keel - 6 kts. 
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Figure 70. Sea state vs. wave angle for 37.5x1 bilge keel - 8 kts. 
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Figure 71. Sea state vs. wave angle for 37.5x1 bilge keel - 10 kts. 
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Figure 72. Sea state vs. wave angle for 37.5x1 bilge keel - 12 kts. 
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Maximum sea state vs. wave angle (for all speeds and 2 degree roll angle) 





Figure 73. Sea state vs. wave angle, no bilge keel, 2 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 4 degree roll angle) 





Figure 74. Sea state vs. wave angle, no bilge keel, 4 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 6 degree roll angle) 
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Figure 75. Sea state vs. wave angle, no bilge keel, 6 degrees roll. 
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Figure 76. Sea state vs. wave angle, 150x3 bilge keel, 2 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 4 degree roll angle) 





Figure 77. Sea state vs. wave angle, 150x3 bilge keel, 4 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 6 degree roll angle) 





Figure 78. Sea state vs. wave angle, 150x3 bilge keel, 6 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 2 degree roll angle) 





Figure 79. Sea state vs. wave angle, 150x2 bilge keel, 2 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 4 degree roll angle) 





Figure 80. Sea state vs. wave angle, 150x2 bilge keel, 4 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 6 degree roll angle) 





Figure 81. Sea state vs. wave angle, 150x2 bilge keel, 6 degrees roll. 
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Figure 82. Sea state vs. wave angle, 150x1 bilge keel, 2 degrees roll. 
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Figure 83. Sea state vs. wave angle, 150x1 bilge keel, 4 degrees roll. 
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Figure 84. Sea state vs. wave angle, 150x1 bilge keel, 6 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 2 degree roll angle) 





Figure 85. Sea state vs. wave angle, 75x3 bilge keel, 2 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 4 degree roll angle) 





Figure 86. Sea state vs. wave angle, 75x3 bilge keel, 4 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 6 degree roll angle) 





Figure 87. Sea state vs. wave angle, 75x3 bilge keel, 6 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 2 degree roll angle) 





Figure 88. Sea state vs. wave angle, 75x2 bilge keel, 2 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 4 degree roll angle) 





Figure 89. Sea state vs. wave angle,75x2 bilge keel, 4 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 6 degree roll angle) 





Figure 90. Sea state vs. wave angle, 75x2 bilge keel, 6 degrees roll. 


[oO 





Maximum sea state vs. wave angle (for all speeds and 2 degree roll angle) 





Figure 91. Sea state vs. wave angle, 75x1 bilge keel, 2 degrees roll. 
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Figure 92. Sea state vs. wave angle, 75x1 bilge keel, 4 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 6 degree roll angle) 





Figure 93. Sea state vs. wave angle, 75x1 bilge keel, 6 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 2 degree roll angle) 





Figure 94. Sea state vs. wave angle, 37.5x3 bilge keel, 2 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 4 degree roll angle) 





Figure 95. Sea state vs. wave angle, 37.5x3 bilge keel, 4 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 6 degree roll angle) 





Figure 96. Sea state vs. wave angle, 37.5x3 bilge keel, 6 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 2 degree roll angle) 





Figure 97. Sea state vs. wave angle, 37.5x2 bilge keel, 2 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 4 degree roll angle) 





Figure 98. Sea state vs. wave angle, 37.5x2 bilge keel, 4 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 6 degree roll angle) 





Figure 99. Sea state vs. wave angle, 37.5x2 bilge keel, 6 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 2 degree roll angle) 





Figure 100. Sea state vs. wave angle, 37.5x1 bilge keel, 2 degrees roll. 
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Maximum sea state vs. wave angle (for all speeds and 4 degree roll angle) 





Figure 101. Sea state vs. wave angle, 37.5x1 bilge keel, 4 degrees roll. 


Maximum sea state vs. wave angle (for all speeds and 6 degree roll angle) 





Figure 102. Sea state vs. wave angle, 37.5x1 bilge keel, 6 degrees roll. 
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B. OPERABILITY INDICES 


Operability Indices, as used in this study, can provide a relative 
measure of the ship’s ability to operate under all sea states and speeds. 
From the polar plot for each speed and a given maximum roll angle, the 
area encircled by the resulting curve could be calculated utilizing the 
trapezoidal method by the use of the trapz function in MATLAB. 
Similarly, the maximum area of the circle can also be calculated if the 
ship is assumed to be able to operate in all sea states without exceeding 
the maximum allowed roll angle. The operability index can then be 
calculated by dividing the calculated area of interest by the maximum 
area. Appendix D provides the program used to calculate the operability 


indices. 


By comparing the operability index related toa particular bilge keel 
size and ship’s speed to that of no bilge keel, a relative measure of the 
effectiveness of that particular bilge keel can be determined. Figure 103 
is a plot of the Operability indices for various ship’s speeds, no bilge keel, 
and 2,4, and 6 degree maximum roll angles. Since one objective of this 
study was to determine the effectiveness of various bilge keels on the roll 
characteristics of the ship, a better measure of the effectiveness was 
found by presenting all bilge keel operability indices with respect to the 


no bilge keel operability indices. In other words, the operability indices 


1g 


of the bilge keel under consideration were divided by the respective 
operability indices for the ship with no bilge keel. If the two values were 
exactly the same, the result would be 1.0. By multiplying the result with 
the corresponding no bilge keel operability index value, one could 
determine the actual respective operability index value. Figures 104 to 


112 are plots of these relative values for the range of ship’s speed. 


Further examination of the operability indices provided some 
useful points. First, the operability indices seem lowest at 2 and 8 knots 
for the no bilge keel ship with the absolute minimum at 8 knots. Six 
knots and speeds above ten knots provide the speeds best suited to allow 
the ship the highest amount of operability. This results in the bilge keels 
having less of an effect at these speeds and the greatest effect at 2 and 8 
knots. As can be seen, any bilge Keel arrangement still significantly 
improves the operability of the ship at most any ship speed. It appears 
that shortening the length of the bilge keel by half while maintaining the 
same width reduces the effectiveness by about 0.2. If the length is 
maintained constant while the width is reduced in size by one foot, the 
overall reduction in effectiveness appears to be approximately 0.2 also. 
These observations are generalities and can vary by a greater or less 


extent based on the different bilge keels and roll angles examined. 
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Figure 103. Operability Index for no bilge keel. 
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Figure 104. Operability Index for 150x3 bilge keel. 
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Figure 105. Operability Index for 150x2 bilge keel. 
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Figure 106. Operability Index for 150x1 bilge keel. 
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Figure 107. Operability Index for 75x3 bilge keel. 
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Figure 108. Operability Index for 75x2 bilge keel. 





Operability Index for 75x1 bilge keel and 2,4, and 6 deg roll angle 
1.35 


ies 


1.25 


lee 





1.15 


1.1 


Operability Index Comparison 


1.05 


0.95 
0 2 4 6 8 10 12 


Speed (kts) 
Figure 109. Operability Index for 75x1 bilge keel. 
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Figure 110. Operability Index for 37.5x3 bilge keel. 
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Figure 111. Operability Index for 37.5x2 bilge keel. 
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Figure 112. Operability Index for 37.5x1 bilge keel. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


AS was seen, not only the size of a wave incident to a ship but also 
the ship’s speed can produce large roll effects. Wave direction also 
played a large part on the magnitude of the ship’s roll. The data provided 
favorable results regarding the use of bilge keels to reduce the 
magnitude of the ship’s roll. These results, though, are based on theory 
and historical results. So, the possibility exists that actual bilge keel 
performance could and probably will vary from the results found in this 


study. The following conclusions can be drawn: 


1. AS was expected, the maximum allowed sea state increased as the 


maximum allowed roll angle was increased. 


2. As ship’s speed was increased, the maximum allowed sea state 


correspondingly decreased. 


3. Beam seas seemed to have a significant effect on the roll 
characteristics of the ship and roll was zero for sea directions relating to 
the bow or stern of the ship. The largest effect, however, on roll angle 


came from quartering seas with quartering seas aft the beam having the 
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largest effect. Thus, aft quartering seas being the limiting position for 


roll angle limits. 


4. Significant roll dampening can be applied by the use of bilge keels. 


5. The operability indices seem lowest at 2 and 8 knots for the no bilge 
keel ship with the absolute minimum at 8 knots. Thus, bilge keels will 


have the greatest effect at these speeds. 


6. Any bilge keel arrangement can significantly improve the operability of 


the ship. 


7. Speeds below six knots and above ten knots provided the ship the 
highest amount of operability with no bilge keel. Bilge keels still proved 


to significantly decrease the roll motion of the ship even at these speeds. 


8. Shortening the length of the bilge keel by half while maintaining the 


same width appeared to reduce the effectiveness by approximately 0.2. 


9. If the bilge keel length is maintained constant while the width is 


reduced in size by one foot, the overall reduction in effectiveness appears 


to be approximately 0.2 also. 
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B. RECOMMENDATIONS 
To continue the research on this topic, the following 


recommendations are made: 


1. This study investigated the effects of bilge keels on the roll 
characteristics of the ship. Several other methods of providing roll 


dampening exist and could also be explored. 


2. Once the placement of the radar masts, or other equipment of 
concern, have been determined, motion characteristics at these locations 
could be determined. These motions would result not only because of 


roll but also because of the other modes of motion of the ship. 


3. Once the maximum allowable roll angle for proper operation of the 
masts and/or equipment of concern has been determined, the minimum 
required bilge keel size can then be determined by the same methods of 


this study. 
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APPENDIX A. SAMPLE SHIPMO DATA INPUT FILE. 


Initial for T-agos 

ae eno 1 ft 0 0 1 O O17 0 
224.00 30.00 3.00 3.00 22.96 5.28 12.00 
224.0000 1.9905 32.1740 0.00000E+00 0.00000E+00 0.0000 
18.7500 -18.7500 1.0000 
1 97.0000 0.0000 0 
0.0000 -15.0000 

6 87.3000 0.0000 0 
0.0000 -15.0000 

0.5208 -12.0000 

1.8750 -9.0000 
3.0458 -6.0000 
5.0000 -3.0000 

6.5000 0.0000 

6 77.6000 0.0000 0 
0.0000 -15.0000 

3.1563 -12.0000 

5.8125 -9.0000 
8.5417 -6.0000 
9.9375 -3.0000 

10.9790 0.0000 

6 67.9000 0.0000 0 
0.0000 -15.0000 

5.4479 -12.0000 
10.0000 -9.0000 
12.4896 -6.0000 
13.6250 -3.0000 

14.7083 0.0000 

6 58.2000 0.0000 0 
0.0000 -15.0000 
7.9792 -12.0000 
14.4375 -9.0000 
15.6146 -6.0000 
16.6875 -3.0000 

17.5417 0.0000 

6 48.5000 0.0000 0 
0.0000 -15.0000 
11.0729 -12.0000 
17.0000 -9.0000 
17.9479 -6.0000 
19.1667 -3.0000 
19.5625 0.0000 

6 38.8000 0.0000 0 
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0.0000 
14.4896 
18.6250 
Lo oon 
20.6250 
20.8750 


-15.0000 


-12.0000 
-9.0000 
-6.0000 
-3.0000 

0.0000 


6 19.4000 0.0000 0 
0.0000 -15.0000 


19.0938 
19°9375 
20.6354 
21.5000 
21.5000 


- 12.0000 
-9.0000 
-6.0000 
-3.0000 

0.0000 


6 9.0000 90.0000 0 


0.0000 
19.5833 
20.0625 
20.8229 
21.5000 
21.5000 


- 15.0000 


-12.0000 
-9.0000 
-6.0000 
-3.0000 

0.0000 


6 -19.4000 0.0000 0 


19.5833 
20.1875 
20.8229 
21.5000 
21.5000 
21.5000 


- 15.0000 
-12.0000 
-9.0000 
-6.0000 
-3.0000 
0.0000 


6 -38.8000 0.0000 0 
0.0000 -15.0000 


19.1354 
20.1250 
20.7500 
21.5000 
21.5000 


-12.0000 
-9.0000 
-6.0000 
-3.0000 

0.0000 


6 -58.2000 0.0000 0 


0.0000 
19.1250 
19.6250 
20.1146 
20.5000 
21.5000 


-11.6250 


-9.0000 
-7.5000 
-6.0000 
-3.0000 

0.0000 


6 -67.9000 0.0000 0 


0.0000 
3.1250 
13.7500 
19.1042 


-9.2500 

-9.0000 
-7.5000 
-6.0000 
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19.9375 -3.0000 
21.4792 0.0000 

6 -77.6000 0.0000 O 
0.0000 -6.8750 

8.1667 -6.0000 

18.2500 -4.5000 
19.2500 -3.0000 
20.3750 -1.5000 
21.3750 0.0000 

6 -87.3000 0.0000 0 
0.0000 -4.7500 

3.0000 -4.6250 

6.0000 -4.2500 
17.2500 -3.0000 
18.8750 -1.5000 
20.6250 0.0000 

6 -97.0000 0.0000 0 
0.0000 -3.0000 

6.0000 -2.7500 
12.0000 -2.5000 
15.0000 -2.3750 
18.0000 -0.6250 
18.9792 0.0000 

1 -103.8759 0.0000 O 
0.0000 -2.0000 
5.2800 0.0000 
9999.0000 0.0000 0.000 
1.0000 0.2000 2.0000 0.1000 13.5025 13.5025 0.0000 
0.0000 180.0000 22.5000 
1.9000 
0.0 
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APPENDIX B. BODY PLANS FOR T-AGOS CLASS 4 SHIP. 
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APPENDIX C. SAMPLE DATA AND OUTPUT FILE. 


% Sets up required data sets for sea state, roll angle, and angle on ship 
% Speed = 0-12 kt., Sea State = 1:8. 150x1 bilge keel 
clear 


th1=[0.0 22.5 45.0 67.5 90.0 112.5 135.0 157.5 180.0]; 
ss=[1234567 8]: 

speed=[0 2 4 6 8 10 12]; 
w=[0.0000,0.1190,0.2590,0.3350,0.1930,0.2490,0.2340,0.1390,0.0000: 
0.0000,0.8550, 1.4100, 1.3900,0.8300, 1.1400, 1.0700,0.7080,0.0000; 
0.0000, 1.5700,2.3600,2.2400, 1.5200, 1.9400, 1.9100, 1.2900,0.0000; 
0.0000,2.2900,3.2400,3.0500,2.3100,2.7500,2.7700, 1.9600,0.0000; 
0.0000,3.7500,4.9600,4.7500,4.0600,4.5000,4.5600,3.4200,0.0000;: 
0.0000,4.5800,6.0000,5.8200,5.1800,5.5700,5.5900,4.2600,0.0000; 
0.0000,5.2600,6.8700,6.7600,6.2800,6.6100,6.5100,4.9200,0.0000: 
0.0000,5.4400, 7.2000, 7.1700,6.7300,7.0300,6.8500,5.1400,0.0000: 
0.0000,0.1760,0.3400,0.3730,0.1900,0.2220,0.1790,0.0991,0.0000: 
0.0000, 1.1700,1.7200, 1.5000,0.8200, 1.0300,0.8610,0.5190,0.0000; 
0.0000,2.0000,2.7400,2.3700, 1.5000, 1.7600,1.5700,0.9640,0.0000;: 
0.0000,2.7400,3.6000,3.1900,2.2800,2.5500,2.3800, 1.5600,0.0000;: 
0.0000,4.1000,5.2000,4.8100,4.0400,4.3400,4.2300,3.0800,0.0000: 
0.0000,4.8200,6.1200,5.8200,5.1300,5.5000,5.3900,4.0000,0.0000: 
0.0000,5.3400,6.9400,6.7000,6.2300,6.5400,6.4100,4.7800,0.0000: 
0.0000,5.5300,7.2400,7.0800,6.6800,7.0300,6.8100,5.0400,0.0000; 
0.0000,0.2900,0.4660,0.4180,0. 1860,0.2000,0. 1420,0.0740,0.0000; 
0.0000, 1.6500,2.1000, 1.6000,0.7990,0.9320,0.7070,0.3970,0.0000; 
0.0000,2.5600,3.1500,2.4800, 1.4400,1.5800, 1.3000,0.7390,0.0000; 
0.0000,3.2700,3.9900,3.2500,2.1900,2.3100,2.0100, 1.2300,0.0000; 
0.0000,4.4400,5.4700,4.8100,3.9200,4.1000,3.8400,2.6900,0.0000; 
0.0000,5.0300,6.3000,5.7600,5.0600,5.3100,5.0800,3.7100,0.0000; 
0.0000,5.4800,6.9700,6.6000,6. 1100,6.4800,6.2200,4.5600,0.0000; 
0.0000,5.6700,7.2300,6.9600,6.5700,6.9600,6.6600,4.8600,0.0000; 
0.0000,0.5900,0.7010,0.4760,0.1840,0.1810,0.1150,0.0572,0.0000; 
0.0000,2.3800,2.6000, 1.7200,0.7860,0.8600,0.6050,0.3190,0.0000; 
0.0000,3.2100,4.5700,2.5500, 1.3700, 1.4400, 1.1000,0.5940,0.0000; 
0.0000,3.8100,4.3000,3.2500,2.0300,2.0800, 1.7000,0.9630,0.0000; 
0.0000,4.7500,5.6300,4.6600,3.6400,3.7700,3.3800,2.2300,0.0000; 
0.0000,5.2400,6.3400,5.5300,4.7700,5.0400,4.7000,3.3000,0.0000; 
0.0000,5.6200,7.0000,6.4000,5.8800,6.2600,5.9300,4.2500,0.0000; 
0.0000,5.7500,7.2300,6.7500,6.3600,6.7900,6.4200,4.6100,0.0000; 
0.0000, 1.5500, 1.2400,0.5460,0.1810,0.1650,0.0959,0.0455,0.0000;: 
0.0000,3.5700,3.4500, 1.9100,0.7850,0.8000,0.5270,0.2640,0.0000;: 
0.0000,4.2700,4.3500,2.7400, 1.3600, 1.3500,0.9770,0.5000,0.0000;: 
0.0000,4.6900,4.9800,3.4000, 1.9900, 1.9600, 1.5200,0.8180,0.0000; 
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0.0000,5.4300,6.0700,4.6800,3.5200,3.5700,3. 1200, 1.9700,0.0000; 
0.0000,5.7700,6.7000,5.4700,4.6100,4.8200,5.6200,3.0200,0.0000; 
0.0000,6.0700, 7.2300,6.2400,5.7300,6. 1400,5.7400,4.0100,0.0000; 
0.0000,6.1700,7.4400,6.5700,6.2100,6.6700,6.2800,4.3900,0.0000; 
0.0000,0.9200,2.5200,0.6370,0.1780,0.1500,0.0811,0.0372,0.0000; 
0.0000,3.1600,4.7300,2.1400,0.7810,0.7440,0.4640,0.2230,0.0000; 
0.0000,3.9800,5.4900,2.9900, 1.3700, 1.2700,0.8690,0.4280,0.0000; 
0.0000,4.4900,6.0100,3.6600,2.0000, 1.8700, 1.3700,0.7070,0.0000; 
0.0000, 5.2800,6.9400,4.8600,3.5300,3.4900,2.9500, 1.8100,0.0000; 
0.0000,5.6700,7.4400,5.5700,4.6000,4.7700,4.2800,2.8700,0.0000; 
0.0000,5.9700,7.8600,6.2800,5.7100,6.1200,5.6500,3.9200,0.0000; 
0.0000,6.0700,7.9800,6.6000,6.2100,6.6700,6.2300,4.3300,0.0000; 
0.0000,0.5980,2.3000,0.7630,0.1750,0.1380,0.0698,0.0310,0.0000; 
0.0000,2.0400,5.2400,2.4300,0.7770,0.6950,0.4140,0.1930,0.0000; 
0.0000,2.7200,6. 1000,3.2900, 1.3600, 1.1900,0.7790,0.3720,0.0000; 
0.0000,3.2100,6.6900,3.9200,2.0000, 1.7700, 1.2400,0.6170,0.0000; 
0.0000,4.0400, 7.5300,5.0600,3.5400,3.4000,2.7800, 1.6600,0.0000; 
0.0000,4.4900,8.0000,5.7700,4.6100,4.7200,4.1700,2.7500,0.0000; 
0.0000,4.8400,8.3900,6.4400,5.7000,6.1200,5.6000,3.8600,0.0000; 
0.0000,4.9600,8.5500,6.7500,6.2000,6.6900,6.2100,4.2900,0.0000]; 


% Utilizes a spline function to smooth data 


for j= 1:7 

cir=T; 

en=(j-1)*8 + 1; 

beg=en+7; 

for i=beg:-1l:en 
t=0:180; 
w1(i,:)=spline(th 1,wi{i,:),t); 
w1p(ctr,:)=w1({i,:); 

ctr=ctr+1; 

end 

figure(j) 

plot(t,w1 p(:,:)) 


xlabel("‘Wave Angle (degrees - 150x1 bilge Keel)’) 

ylabel('Roll Angle (degrees)’) 

s=speed(j); 

title(["Roll angle vs. wave angle (',int2str(s),' kts., sea states (ss) 1 to 
8)'}) 

legend('ss8','ss7','ss6','ss5','ss4','ss3','ss2','ss1',-1) 
end 
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% Calculates maximum roll angle (2,4,6 degrees) for angle on ship. 
ctr1=8; 
etrZ2=1; 
for jl=1:7 
beg=(j1-1)*8 + 1; 
en=beg+7; 
ctr=1; 
for k=beg:en 
w la(ctr,:)=w1(k,:); 
ctr=ctr+1; 
end 


ra2=8.0*ones(1,181); 
ra4=8.0*ones(1,181); 
ra6=8.0*ones(1, 181); 
[rl,cl]=size(wla); 
ssn=1.00:0.01:8.00; 
ssnc=size(ssn,2); 


for i=I:cl 

wne=spline(ss,w1la(:,i),ssn); 

if wn(1,1) > 2.0 
ra2(i)=0.0; 

end 

if wn(1,1) > 4.0 
ra4(i)=0.0; 

end 

if wn(1,1) > 6.0 
ra6(ij=0.0; 

end 


for j=1:(ssnc- 1) 


if wn(1,j) <= 2.0 
k=j+1; 
if wn(1,k) > 2.0 
ra2(i)=ssn(j); 
end 
end 


if wn(1,j) <= 4.0 
k=j+1; 
if wn(1,k) > 4.0 
ra4(ij=ssn{j); 
end 
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end 


if wn(1,j) <= 6.0 
k=j+1; 
if wn(1,k) > 6.0 
ra6(i)=ssn{j); 
end 
end 
end 
end 


% Sets up roll angle and theta matrices for plotting 


for i=1:181 
Z=i+ 180; 
th3(iJ=t(i)*pi/ 180; 
th3(z)=t(i)*pi/180 + pi; 
j=182-i; 
ra2(z)=ra2(j); 
ra4(z)=ra4(j); 
ra6(z)=ra6(j); 

end 


% Prints out polar plot for sea state limits for roll angles of 2,4, and 
% 6 degrees for each speed. 

figure(ctr1) 

polar(th3,ra6,'r') 

hold on 

polar(th3,ra4,'b’) 

polar(th3,ra2,'g’) 

legend('g','2 degrees','b','4 degrees’,'r','6 degrees’,-1) 

xlabel('150x1 bilge keel’) 

hold off 

s2 = speed({j1); 

title([/Maximum sea state vs. wave angle (for ',int2str(s2),’ kt. and 2,4,6 
degree roll angle)']) 

pra2(ctr2,:)=ra2(1,:); 

pra4(ctr2,:)=ra4(1,:); 

pra6(ctr2,:)=ra6(1,:); 

ctrl=ctr1+1; 

ctr2=ctr2+1; 

end 


%Prints out polar plots for sea state limits for all speeds for 2,4, and 
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%6 degrees roll angles. 

figure(15) 

polar(th3, pra2(1,:),'y’) 

legend('0 kts',- 1) 

hold on 

polar(th3,pra2(2,:),'m’) 

polar(th3,pra2(3,:),'c’) 

polar(th3,pra2(4,:),'r) 

polar(th3, pra2(5,:),'g') 

polar(th3,pra2(6,:),'b’) 

polar(th3, pra2(7,:),'w’) 

Xlabel('150x1 bilge keel’) 

title("Maximum sea state vs. wave angle (for all speeds and 2 degree roll 
angle)') 

legend(‘y’,'0 kts','m','2 kts','c',’4 kts','b',/10 kts','w’,'12 kts',-1) 
hold off 


figure(16) 

polar(th3,pra4(1,:),'y’) 

hold on 

polar(th3, pra4(Z, :),'m’') 

polar(th3,pra4(3,:),'c’) 

polar(th3,pra4(4,:),'r’) 

polar(th3,pra4(5,:),'g') 

polar(th3, pra4(6,:),'b') 

polar(th3, pra4(7,:),'w’) 

xlabel('150x1 bilge keel’) 

title(Maximum sea state vs. wave angle (for all speeds and 4 degree roll 
angle)') 

legend(‘y','O kts','m','2 kts’,'c','4 kts’,'b’,/10 kts','w','12 kts’',-1) 
hold off 


figure(17) 

polar(th3,pra6(1,:),'y’) 

hold on 

polar(th3,pra6(2,:),'m') 

polar(th3, pra6(3,:),'c’) 

polar(th3, pra6(4,:),'r’) 

polar(th3, pra6(5,:),'g’) 

polar(th3, pra6(6,:),'b’) 

polar(th3, pra6(7,:),'w’) 

hold off 

xlabel('150x1 bilge keel’) 

title(Maximum sea state vs. wave angle (for all speeds and 6 degree roll 
angle)’ 

legend(‘y','0 kts','m’','2 kts','c’,'4 kts','b’,'10 Kts’,'w’','12 kts',-1) 
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figure(1) 

print -depsc2 g150x10.eps 
figure(2) 

print -depsc2 g150x12.eps 
figure(3) 

print -depsc2 g150x14.eps 
figure(4) 

print -depsc2 g150x16.eps 
figure(5) 

print -depsc2 g150x18.eps 
figure(6) 

print -depsc2 g150x110.eps 
figure(7) 

print -depsc2 g150x112.eps 
figure(8) 

print -depsc2 p150x10.eps 
figure(9) 

print -depsc2 p150x12.eps 
figure(10) 

print -depsc2 p150x14.eps 
figure(1 1) 

print -depsc2 p150x16.eps 
figure(12) 

print -depsc2 p150x18.eps 
figure(13) 

print -depsc2 p150x110.eps 
figure(14) 

print -depsc2 p150x112.eps 
figure(15) 

print -depsc2 p150x12d.eps 
figure(16) 

print -depsc2 p150x14d.eps 
figure(17) 

print -depsc2 p150x16d.eps 


% Calculates areas under 2,4, and 6 degree curves to find operability 
index 
for i=1:7 
a15012(i)=trapz(th3,pra2(i,:)); 
a15014(i)=trapz(th3,pra4ii,:)); 
a15016(ij)/=trapz(th3, praG(i,:)); 
end 
diary on 
al5012 
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al5014 
al5016 
diary off 


105 


106 





APPENDIX D. OUTPUT FILE FOR OPERABILITY INDICES. 


% Operability index for no and all bilge keel sizes 


sp=[0 2 4 68 10 12]: 


oi=[20.5352 19.6437 20.9513 21.8156 18.7591 22.2139 


23.1490; 


27.5689 
30.2916 
40.6990 
50.2655 
50.2655 
31.5754 
90.1601 
90.2655 
23.6220 
33.8953 
43.4954 
33.1574 
50.2655 
50.2655 
27.0449 
41.0418 
50.2655 
22.0652 
30.9436 
36.6973 
28.3462 
43.5711 
50.2655 
24.0074 
34.6294 
45.0302 
21.3010 
29.3208 
33.6042 


25.4916 
28.4353 
40.6389 
50.2655 
50.2655 
31.4809 
50.1964 
50.2655 
23.6932 
33.9282 
43.5624 
33.1825 
50.2655 
50.2655 
26.8404 
40.9971 
50.2655 
22.1496 
30.9119 
36.9689 
28.0363 
43.6137 
50.2655 
24.0688 
34.7334 
45.0246 
212152 
28.9641 
33.5519 


26.3032 
29.7059 
40.4951 
50.2655 
50.2655 
31.0794 
49.9981 
50.2655 
23.5707 
34.0430 
43.7010 
33.0098 
50.2655 
50.2655 
26.9161 
40.9294 
50.2655 
22.1688 
30.7803 
37.3480 
27.6600 
43.4905 
50.2655 
23.9180 
34.8465 
45.2023 
21.3695 
28.7246 
33.9582 


27.0467 
30.6955 
40.1558 
50.2655 
50.2655 
30.4644 
49.5953 
90.2655 
23.6384 
33.5065 
44.2738 
32.4087 
50.2655 
50.2655 
26.0630 
41.1273 
50.2655 
22.0700 
30.6469 
38.3693 
27.1524 
43.4775 
50.2655 
23.6660 
35.1848 
45.7350 
21.3614 
28.6705 
34.9195 


24.5390 
28.4974 
39.2916 
50.2655 
50.2655 
29.0074 
49.2371 
50.2655 
22.6258 
32.3755 
43.9998 
30.8072 
50.2655 
50.2655 
26.1045 
41.0348 
50.2655 
21.2421 
29.1149 
37.8269 
26.0138 
42.7571 
50.2655 
22.7242 
oi o 
45.5974 
20.2741 
27.1758 
33.6370 


% changes operability index to % for no bilge keel 


max=50.2665; 


for i=1:3 


pooi(i,:)=oi(i,:). /max; 
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27.0156 
30.8023 
36.7968 
50.2655 
50.2655 
28.0754 
47.6217 
90.2655 
Zoe oe 
32.3211 
43.0580 
ZFS Oe 
49.3792 
50.2655 
24.8517 
38.2164 
90.2655 
21.4560 
29.6535 
37.2617 
25.6165 
41.5179 
50.2655 
22.4861 
33.0551 
44.6895 
21.2836 
28.4457 
33.9735 


Zo. 117: 
32.1364; 
$5.22 76: 
50.2655; 
50.2655; 
28.4607; 
46.2013; 
50.2655; 
24.0342; 
Oo ODO. 
42.8258; 
29.5969; 
47.5442; 
50.2655; 
2587132; 
38.5906; 
49.1335; 
23.5047; 
ol. ia: 
37.9979; 
26.4812; 
40.6086; 
50.1873; 
2490 154 
33.9994; 
44.3698; 
Zo 1232. 
30.0060; 
35.0874]; 


end 
% compares operability index of no bilge Keel to different bilge keels 


for i=4:3:30 
=i+1; 
m=i+2; 
for j=1:7 
pooi(i,j)=oi(i,j). /oi(1,j); 
pooi(l,j)=oi(1,j). /oi(2,j); 
pooi(m,j)=oi(m,j). /oi(3,j); 
end 
end 


%Printing out results 


figure(1) 

plot(sp, pooi(1,:),sp,pooi(2,:),sp,pooi(3,:)) 

xlabel(‘Speed (Kts)’) 

ylabel(‘Operability Index’) 

title('Operability Index for no bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg’,'4 deg’,'6 deg’,-1) 

print -depsc2 oinbk.eps 

figure(2) 

plot(sp,pooi(4,:),sp,pooi(5,:),sp,pooi(6,:)) 

xlabel('Speed (Kts)’) 

ylabel(‘Operability Index comparison’) 

title(‘(Operability Index for 150x3 bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg’,'4 deg','6 deg’,- 1) 

print -depsc2 0i1150x3.eps 

figure(3) 

plot(sp, pooi(7,:),sp,pooi(8,:),sp,pooi(9,:)) 

xlabel('‘Speed (Kts)') 

ylabel(‘Operability Index Comparison’) 

title((Operability Index for 150x2 bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg','4 deg','6 deg’,-1) 

print -depsc2 0i1150x2.eps 

figure(4) 

plot(sp, pooi(10,:),sp,pooi(11,:),sp,pooi(12,:)) 

xlabel(‘Speed (kts)') 

ylabel(‘Operability Index Comparison’) 

title((Operability Index for 150x1 bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg’,'4 deg’,'6 deg’,-1) 

print -depsc2 0i1150x1l.eps 

figure(5) 
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plot(sp,pooi(13,:),sp,pooi(14,:),sp,pooi(15,:)) 

xlabel(‘Speed (kts)’) 

ylabel(‘Operability Index’) 

title((Operability Index for 75x3 bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg’,'4 deg’,'6 deg’,-1) 

print -depsc2 0175x3.eps 

figure(6) 

plot(sp, pooi(16,:),sp,pooi(17,:),sp,pooi(18,:)) 

xlabel(‘Speed (kts)’) 

ylabel(‘Operability Index Comparison’) 

title(Operability Index for 75x2 bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg’,’4 deg’,’6 deg’,-1) 

print -depsc2 0i175x2.eps 

figure(7) 

plot(sp, pooi(19,:),sp,pooi(20,:),sp, pooi(2 1,:)) 

xlabel('Speed (Kkts)’) 

ylabel(‘'Operability Index Comparison’) 

title(Operability Index for 75x1 bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg’,'4 deg’,'6 deg’,-1) 

print -depsc2 oi75x1.eps 

figure(8) 

plot(sp, pooi(22,:),sp,pooi(23,:),sp,pooi(24,:)) 

xlabel(‘Speed (Kts)’) 

ylabel(‘Operability Index Comparison’) 

title('Operability Index for 37.5x3 bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg’,'4 deg’,'6 deg’,-1) 

print -depsc2 01375x3.eps 

figure(9) 

plot(sp, pooi(25,:),sp,pooi(26,:),sp,pooi(27,:)) 

xlabel('Speed (Kts)’) 

ylabel(‘'Operability Index Comparison’) 

title((Operability Index for 37.5x2 bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg’,'4 deg’,’6 deg’,-1) 

print -depsc2 01375x2.eps 

figure(1 0) 

plot(sp,pooi(28,:),sp,pooi(29,:),sp, pooi(30,:)) 

xlabel('Speed (kts)’) 

ylabel(‘Operability Index Comparison’) 

title((Operability Index for 37.5x1 bilge keel and 2,4, and 6 deg roll angle’) 
legend('2 deg’,'4 deg’,'6 deg’,-1) 

print -depsc2 01375x1.eps 


HOS 
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